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Abstract. The problem of two fixed centers is a classical integrable problem, 
stated and integrated by Euler in 1760. The integrability is due to the unex- 
pected first integral G. Some straightforward generalizations of the problem still 
have the generalization of G as a first integral, but do not possess the energy 
integral. We present some numerical integrations suggesting that in the domain 
of bounded orbits the behavior of these a priori non hamiltonian systems is very 
similar to the behavior of usual quasi-integrable systems. 

The equations. Euler 's problem in the plane (see Figure 1) is defined by the 
system of differential equations 

X ^ -a{xA,y)xA - b{xB,y)xB, i) = -a{xA,y)y - b{xB,y)y. (1) 

The two fixed centers are the points (1,0) and (—1,0), and the moving particle 
is the point (x, y). We have set xa = x — 1, xb — x + 1, 

a{^,v)=niA{e+v'r^^^ b{Lv)^mB{e+vT^^^- (2) 

The problem can be defined in the 3-dimensional space in the same way, and is 
also integrable, as was noticed by Euler. However, we will restrict ourselves to 
the planar case. 

The first step in Eulcr's integration was to exhibit two independent first 
integrals of the motion. One is the energy 

H={x^+ f)/2 - h{xA, y) - b{xB, y), 

with a(f,r;) = mA(C^ + ?7^)"^/^, Kt.v) = "^b(6^ + V^T^^'^- We will call the 
second one Euler's integral: 

G = GaGb ~ 2a{xA,y)xA + 2b{xB,y)xB, 
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with Ca = xaV — yx, Cb = xbV — yx. Eulcr continued the integration, ehmi- 
nating the second derivatives in (1) using the first integrals, and separating the 
variables. 

Our generalization is simply to consider system (1) in the case where a 
and b are any homogeneous functions of degree —3. Indeed, we want to put a 
little restriction on these homogeneous functions. We will suppose that both 
differential forms £,a{£,,'r]){Cdri — rjdS^) and ria{£^,r]){^dri — r/rf^) arc exact forms 
on the plane minus the origin, and that the same is true when we change a 
in b. This hypothesis comes from the study of the problem of one fixed center 
(see [2] and [3]). It is not a strong restriction: the forms are already closed, so 
the condition on each function is just the cancellation of two scalar quantities, 
namely the integrals of both forms on a closed path around the origin. 

It can be shown that any function 0(^,77) satisfying the above conditions 
comes from a function A{^,ri) homogeneous of degree 1 as follows. Let us 
denote by A^, Ar^ the first derivatives of the function A and by A^^, A^,,, 
the second derivatives. Then 

a{^,ri) = V~^A^^ = -rS"'^i^ = r^Ann- 

The function a(^, ry) in Euler's case (2) is obtained in this way from the function 
A = m/i(^^ + jfY/"^. Wc also have = ^a. It was discovered by the first 
author that G persists in the form 

G = GaGb - 2A^{xA,y) + 2B^{xB,y), 

where -B(^, rj) and B^{^, rj) are associated to the function b in the same way as A 
and are associated to a. In general, no integral takes the place of the energy 
integral. 

Quasi-integrability. We report our numerical exploration of these generalized 
Euler's problems, showing three examples that seem to us significant. In all 

cases we met, the result is either escape or quasi-intcgrable behavior. The third 
experiment displays some islands suggesting non integrability. Magnifying the 
neighborhood of a saddle point, a domain of irregular dynamics can be observed. 

The obvious choice for a Poincarc section is to fix the integral G and take, 
for example, j/ = (y > 0). In each case, we show the iterates of some points 
of this Poincare mapping, the central orbit in the section and some typical 
quasipcriodic orbit. All the orbits in a given field of forces have the same value 
of Euler's integral. Since the examples have very large orbits, we have taken 
throughout the numerical experiments a somewhat arbitrary cut-off criterion 
given by the value of any coordinate or velocity greater than a thousand. In all 
the examples, we have taken A(^, rj) = (5^^ —^r]+5r]'^)/10r, where r = + rf, 
which corresponds to a(^, ry) = (5^^ + 3^7/ + 5?7^)/10r. 

First example. Figures 2 to 4 correspond to B{£^,rj) = — (^'^ — 3^ry^)/4r^ + r, 
and thus rf) = —2(3^77^ — ^^)/r^ + . The Poincare map is close to a linear 
map, on a whole domain delimited by the escape criterion. This is quite strange. 
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An explanation for this phenomena conies from geometrical considerations. We 
have chosen the plane as the domain for the motion, but there is a natural 
bigger domain for this kind of systems. It is the manifold of half lines drawn 
from the origin in a 3— dimensional vector space. Our plane is from this point 
of view just one half of the natural domain, a hemisphere chosen arbitrarily. 
Escaping orbits appear as orbits cut by the boundary of the hemisphere (in the 
classical Kcplcr problem, hyperbolas appear in the same way as cut ellipses). 
The theoretical grounds for this remark may be found in [4] . 

Second example. In Figures 5 to 7, we have chosen B = 4(^4 + rj^Y^'^, 
b = 12^^?7^(^^ + ?7^)^^/^. Here the section displays a wide domain with strong 
torsion but we are still very close from an integrable system. This rises the 
question: what are the integrable systems nearby? We know very few cases 
where our generalized Euler problem is integrable, namely the classical case 
and its projective transformations defined in [4] (which correspond for example 
to replace + 77^ in Eq. (2) for b by any homogeneous quadratic expression in 
^, T], and leave a as it is.) Because we needed to get sufficiently many bounded 
orbits we were probably forced to stay close from integrable cases. 

Third example and final comments. In Figures 8 to 11, we have chosen 

B = (3^^ — £,V + 3ri'^)/3r, b = (C^ + + rf')/r^ . Here the system behaves 
as a typical conservative system close to an integrable system. We are more 
accustomed to observe this in the class of hamiltonian systems, and one can 
argue that maybe the system is hamiltonian for some symplectic form. We do 
not believe so, and rather relate this quasi-integrability to KAM theory applied 
to reversible systems (see [5], Theorem 2.9). Our systems are clearly reversible. 
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Figure 1: In Euler problem as well as in the generalizations we present, a particle 
at the position {x, y) evolves under the action of two centers A and B, with 
respective coordinates (1,0) and (—1,0). 
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Figure 2: Example 1; B = -(^^ _ 3^^2-)/4^2 ^ j.. q ^ 5.50433086. Poincare 
section {y = 0,y > 0), from the central periodic orbit to the neighborhood of 
the last torus before cut-off. 
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Figure 3: Example 1 again. Configuration space of the periodic orbit: x = 
2.352375, x = 0.3621675 and y = 1.050626. 
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Figure 4: Example 1. The configuration space of the last orbit shown in the 
Poincarc section just before our escape criterion is satisfied: x = 3.469875, 
X = 0.3621675 and y = 0.99058834. 




Figure 5: Example 2; B = -{S,^-^rf )/'ir'^+r; G = 10.9200094. Same Poincare 
section as Figure 2, showing again a phase space foliated by invariant tori up to 
the escape orbit. 
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Figure 6: Example 2. The configuration space of the periodic orbit with initial 
conditions x = 2.6, x = 0.2934 and y = 0.8249589. 
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Figure 7: Example 2. Top: A typical torus at a; = 6.1, x = 0.2934 and 
y = 0.3290253; Bottom: A detail of the last torus before escape with initial 
conditions x = 13.6, x = 0.2934 and y = 0.145976133. 
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Figure 8: Example 3; B = (3^^ -^r? + 3r/2)/3r; G = 1.55230255. Top: Poincare 

section as in Figure 2 starting at the central periodic orbit up to escaping 
orbits; Bottom: The configuration space of the central periodic orbit x = 4.0005, 
X = 0.191860239 and y = 0.191860239. 
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Figure 9: A magnification of Example 3 Poincare Section showing four islands 
of a chain with fourteen islands. 
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Figure 10: A torus very near the separatrix and its section for Example 3. 
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Figure 11: Example 3. Top: The neighborhood of the islands of the Poincare 
section; Bottom: Its stochastic zones indicating the absence of integrability. 
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